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DESIGN OF ROUND-TRIP TRAJECTORIES TO NEAR-EARTH 
ASTEROIDS UTILIZING A LUNAR FLYBY 

Sonia Hernandez* and Brent W. Barbee 1 


There are currently over 7, 700 known Near-Earth Asteroids (NEAs), and more 
are being discovered on a continual basis. Current models predict that the actual 
order of magnitude of the NEA population may range from 10 5 to 10 6 . The close 
proximity of NEA orbits to Earth's orbit makes it possible to design short duration 
round-trip trajectories to NEAs under the proper conditions. In previous work, 59 
potentially accessible NEAs were identified for missions that depart Earth between 
the years 2016 and 2050 and have round-trip flight times of a year or less. We 
now present a new method for designing round-trip trajectories to NEAs in which 
the Moon’s gravity aids the outbound trajectory via a lunar flyby. In some cases 
this gravity assist can reduce the overall spacecraft propellant required for the 
mission, which in turn can allow NEAs to be reached which would otherwise be 
inaccessible to a given mission architecture. Results are presented for a specific 
case study on NEA 2003 LN 6 . 


INTRODUCTION 

There are currently over 7, 700 known Near Earth Asteroids (NEAs) and more are being discov- 
ered on a continual basis*. NEAs are small (often rocky but occasionally metallic) celestial bodies 
with mean diameters ranging from several meters to several kilometers. Their solar orbits are such 
that they often closely approach Earth’s orbit, and this makes them a potential threat to life on Earth 
by virtue of collision hazards. 1 However, this also presents a unique opportunity for exploration. 
Because NEAs have orbits in close proximity to Earth’s orbit, it is possible to design short duration 
round-trip trajectories to NEAs, which may permit these objects to become destinations for human 
missions. 23 

There are a number of compelling reasons to send missions to NEAs. Asteroids have remained 
largely unchanged in composition and studying them could provide insight into the origins of our 
solar system. Moreover, because an NEA could be on a collision course with Earth in the near 
future, understanding their characteristics is vital. Additionally, a human mission to an NEA would 
be the most ambitious journey ever undertaken by our space program and would renew passion 
for space exploration among both students and the general public. Indeed, the current presidential 
administration has proposed, as part of the Flexible Path program, that such a mission be launched 
during the mid-2020s^. 
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In previous work, a study was performed to identify NEAs that may be accessible for round-trip 
missions that employ a single heavy-lift launch, depart Earth between the years 2016 and 2050, 
and have a total round-trip flight time of less than a year. 4 The algorithm developed for that study 
identified 59 accessible NEAs and 10 marginally inaccessible NEAs. This parametric study was 
performed with two-body dynamics and patched conics between the Earth and the Sun for the 
spacecraft trajectory segments, though full high-precision ephemerides were used for the Earth and 
the NEAs. 

In this paper we present a newly developed method for designing round-trip trajectories to NEAs 
in which the Moon’s gravity aids the outbound segment of the trajectory via a lunar flyby. The goal 
is to perform the flyby in a manner that reduces the overall required propellant mass for the mission, 
which can cause previously inaccessible NEAs to become accessible (assuming no other changes 
to the mission architecture). We demonstrate the efficacy of this methodology by extending the 
work in Reference 4 to include the Moon’s gravity in a three-body model. This is combined with 
an automated algorithm that performs a parametric scan to target round-trip trajectories to NEAs 
utilizing a lunar flyby in a manner that serves to reduce the required propellant mass for the mission. 

This paper is organized as follows. The first section introduces a three-body model, centered at 
Earth and perturbed by the Moon, which serves as a first estimate to design lunar flybys. In the 
next section, we describe the automated algorithm used to design round-trip trajectories to NEAs. 
This algorithm is implemented by combining the model described in the first section with patched 
conics between the Earth and the Sun. The last section deals with implementing Moon’s accurate 
ephemerides into the algorithm and an example trajectory misson is shown for NEA 2003 LN6- 

Disclaimer Regarding NEA Accessibility 

The discussion we present in this paper often hinges on whether a NEA is classified as “ac- 
cessible” or “inaccessible” or “marginally inaccessible” by our algorithms. It is important to note 
that the meaning of “accessible” in the context of human missions to NEAs is currently an area 
of vigorous research and we are making no claims whatsoever regarding whether any NEAs are 
truly accessible for human exploration. Rather, we are attempting to identify potentially accessible 
NEAs from a dynamics point of view, and this what we mean when talk about “accessibility.” How- 
ever, while the trajectory design results may make a particular NEA appear accessible for a given 
mission architecture (mass of crew vehicle, number of launches, type of propulsion system, etc.), 
there are a host of other constraints and considerations (e.g., various operational constraints and 
human health and safety constraints) that must be analyzed in great detail before the accessibility of 
a NEA can be truly and properly evaluated. Those other aspects of human space flight to NEAs are 
quite complex and are thus well outside the scope of this paper. So, to prevent our verbiage from 
becoming unwieldy we will simply talk about NEA “accessibility” in the context of our current 
dynamical analysis without repeating these caveats, under the assumption that the reader has read 
this disclaimer and understands what we mean by “accessible.” 

LUNAR FLYBY IN A THREE-BODY MODEL 

The purpose of this research is to design trajectories to NEAs which utilize a lunar flyby in the 
outbound (Earth departure) segment of the trajectory. The first step in designing the trajectory 
consists of targeting a lunar flyby that will cause the trajectoiy to intercept the NEA. To accomplish 
this, we introduce a dynamics model for the spacecraft which includes the effects of lunar gravity 
in a non-rotating Earth-centered frame. 
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The acceleration experienced by the spacecraft, with Earth as the primary central body and the 
Moon as a perturbing body is 5 
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where r 6 R 3 is the position of the spacecraft with respect to the Earth and r f nl (t) e R 3 is the 
position of the Moon with respect to the Earth at any given epoch. The gravitational parameters of 
the Earth and the Moon are, /. i e and fi m , respectively*. A numerical integration algorithm is used to 
integrate Eq. (1) by defining a state vector 


x = 
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€ R 6 


for the spacecraft with initial conditions xo = [ro vo] T , where vo € R 3 is the velocity at time to. 

We obtain an initial trajectory estimate by modeling the motion of the Moon about the Earth with 
a simplified planar, circular model in which the mean motion of the Moon is 


n m — 

with a m = 384, 400 km (the mean distance between the Earth and the Moon). In this model, the 
position of the Moon, r^(£), at any epoch t is 



r m(f) = “m 


cos(a<) 
sin (at) 
0 


(2) 


The angle a t is measured from the positive x-axis of the Earth-Moon line and a to represents the 
initial location of the Moon when the spacecraft departs Earth (see Figure 1 ). at at any epoch is 


a t = a to + n m t (3) 

For the purposes of this study, we make the assumption that the spacecraft departs from a circular 
Low Earth Orbit (LEO) at an altitude, Hleo^ of 185 km, and travels to the Moon on a Hohmann 
transfer. However, the design of the algorithm is such that any Iileo or an y type of transfer may be 
used. The initial impulse required to arrive at the Moon on this trajectory is 

A VlEO = Vtrans ~ 

= \j^ e {jZEo ~ 25 ^) ” ^ 

= 3.1352 km/s 

We also need to determine where the Moon should be initially located so that the spacecraft 
will fly past the Moon in a manner that provides additional velocity change sufficient to cause the 

*The values for the gravitational parameters are = 3.986 • 10 5 km 3 s~ 2 and = 4.903 • 10 3 km 3 s~ 2 . 
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Figure 1. Schematic of Lunar Flyby in a Model Centered at the Earth and Perturbed by the Moon. 


spacecraft to reach Earth’s Sphere of Influence ( SOI e ). After imparting on the spacecraft located in 
LEO the initial impulse A Vleo> its position and velocity are 


ro 


^LEO COS (0 O ) 

vleo sin(<9 0 ) 


and v 0 = 


— (t; 0 + A Vleo) sin(0 o ) 
(^o + A Vleo) cos (0 o ) 


(5) 


where the angle 9o defines the initial location of the spacecraft from the positive x-axis of the Earth- 
centered inertial frame at time to and tleo = r e + h leo (with r e = 6, 378.137 km). 

An algorithm was created using the model described in Eq. ( 1 ) to search for all the possible angles 
a to such that the lunar gravity assist would cause the spacecraft to escape Earth and reach SOI e . 
The result of this search (for trajectories with Earth departure characteristics as defined in Eq. (4)) 
is that the location of the Moon (at the time at which the spacecraft departs Earth) must be within 
the interval 

6o + 110° < at 0 < 0 0 + 125° (6) 

This means that for initial angles at 0 which are smaller than 110° or greater than 125°, the Moon 
is too far from the spacecraft when it reaches the distance of a m to provide the spacecraft with the 
necessary gravity assist to escape Earth. As can be observed from Figure 2, the spacecraft will 
reach SOI e with different flight paths and velocities depending on the particular value of at 0 within 
the interval given in Eq. (6). Thus, there are families of post-flyby trajectories that may or may 
not reach the NEA even though they do reach SOI e \ merely reaching SOI e is a necessary but 
insufficient condition for subsequent interception of the NEA. 

Note that in the approximated planar, circular model described in this section, the angle of depar- 
ture 6q can vary from 0° < 6o < 360° and a to is varied as given by Eq. (6). When searching for 
trajectories to intercept NEAs, it is necessary first to scan through all the possible combined initial 
conditions of 6q and a to that will produce a trajectory which reaches SOI e . This requires repeat- 
edly numerically integrating Eq. (1), which can be computationally expensive. However, since in 
this system the Moon rotates at a constant rate about the Earth, it is possible to only integrate for 
6o = 0 and any a to as in Eq. (6), which then allows a simple rotation through angle 9o to trans- 
form the trajectory to the desired location. If the state obtained with 9q = 0 is x° = [r° v°] T , the 
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Velocity at SOI of Earth vs. Oq 




Figure 2. (a) Possible Flyby Trajectories that Reach SOI e (b) Velocity Magnitude of 
the Spacecraft Upon Reaching SOI e . 


associated state at any other 6q will be 
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where the transformation matrix R is 


R = 


cos(#o) sin(#o) 0 
— sin(#o) cos(#o) 0 
0 0 1 


(7) 


PATCHED CONICS METHODOLOGY 

In Reference 4, a study was performed to identify NEAs that may be accessible for round-trip 
missions which employ a single heavy-lift launch, depart Earth between the years 2016 and 2050, 
and have a total round-trip flight time of a year or less. The algorithm developed for that study 
identified 59 accessible NEAs and 10 marginally inaccessible NEAs. The parametric analysis per- 
formed in that study utilized two-body dynamics and patched conics between the Earth and the Sun 
for the spacecraft trajectory segments, though full high-precision ephe me rides were used for the 
Earth and the NEAs. 

In this section we extend the work in Reference 4 by designing round-trip trajectories to NEAs in 
which a lunar gravity assist aids the outbound segment of the trajectory. The goal is to perform the 
lunar flyby in a manner that reduces the overall required propellant mass for the mission, which can 
cause NEAs found to be inaccessible in the previous study to become accessible (under the mission 
architecture assumptions of the previous study). 

The profile of a human mission to an NEA utilizing a lunar flyby is shown in Figure 3. The 
round-trip trajectory consists of four stages: Earth departure to lunar flyby, post-flyby trajectory to 
NEA arrival, stay time on the NEA’s orbit, and NEA departure to Earth return. For the purposes of 
this study, the Earth departure date, flight times for each stage, and AV magnitudes were subject to 
the same constraints specified in Reference 4. Each of the four stages are now detailed in turn. 
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Stage 1 Shown in Figure 3(a), this stage of the trajectory is designed in an Earth-centered, inertial 
frame such that the trajectory departs Earth to perform a flyby at the Moon. 

Specifying the date of departure td ePe , the initial departure angle 60 , and the location of the 
Moon at departure a to (see Figure 1) fully defines the initial conditions given in Eq. (5), which 
allows for the integration of Eq. (1). The trajectory is propagated until it reaches SOI e , over 
a time span of tof so i. 

The A Vleo in Eq. (4) corresponds to a characteristic energy for Earth departure, C3 = 
—2.0390 km 2 /s 2 . This value is defined as the square of the hyperbolic excess velocity, 1^, 
with respect to Earth on a hypothetical outbound hyperbola, and is computed from the energy 
equation of an orbit as 6 

E = = ^ - — = ^ =>C 3 = -2.0390 km 2 /s 2 . (8) 

2cim 2^ V oo 2 


Stage 2 When the trajectory reaches SOI e after the flyby, the Earth-centered state x e soi = [r e sai 
is transformed into a Sun-centered frame as shown in Figure 3(b). 
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where x* = [r* Vg] T is the state of the Earth with respect to the Sun at the time when the 
trajectory reaches SOI e , that is, t so i = tde Pe + tof so i, which are computed from precise 
ephemerides files. To simplify the notation, we drop the superscript s from the remaining 
states since they are all in a Sun-centered frame. Note that for this trajectory we use a sim- 
plified planar, circular Earth-Moon model, and therefore x e soi E R 2 , however, x* E R' 3 . We 
have also assumed that the Moon’s trajectory lies in the ecliptic plane. These simplifying 
assumptions are absent in the full version of the algorithm, which is detailed in a subsequent 
section. The true ephemeris of the Moon is utilized in the full version of the algorithm. 

A trajectory from the SOI e crossing to the NEA is now designed via Lambert targeting by 
choosing a time to reach the asteroid from SOI e , tof arra . The trajectory is designed in such 
a way that the post-flyby maneuver, A V so i is minimized. Ideally, A V SO i = 0, a constraint 
that will be added in future work when we perform optimization of the trajectory. This con- 
straint may not always be achievable, but driving AV so i as close to zero as possible is clearly 
beneficial. 

A Vsoi = ||v<iep SO i V so i ||, (10) 

where v^ epsoi is the required velocity at SOI e to intercept the NEA and v so i is the velocity 
that the spacecraft has after the lunar flyby when reaching SOI e . 

At the end of this stage, the spacecraft performs a maneuver, AV arra , to match the orbit of 
the NEA upon arrival. 

AV ar r a ~ II ^flarr ct 1 1 5 

where v aarr is the heliocentric velocity vector of the NEA at the epoch of arrival, which is 
computed from the JPL HORIZONS ephemeris file and v arra is the velocity vector of the 
spacecraft upon reaching the NEA. 


Stage 3 The spacecraft remains with the NEA on its orbit for a selected time span, tof sta y 
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Stage 4 After the stay time at the NEA has elapsed, a trajectory is designed to depart the NEA with 
A Vd epa , and return the spacecraft to Earth with the desired flight time, tof arre . 

A Vdep a = II V<fep a - Va dep ||, (12) 

where v adep is the heliocentric velocity vector of the NEA at the epoch of departure from it, 
which is again computed from the JPL HORIZONS ephemeris file and Vdep a is the velocity 
vector of the spacecraft before departure from it. 

The algorithm also allows for a possible re-entry impulsive maneuver, AV arre if the spacecraft 
exceeds the maximum value v a rr emax = 12 km/ s. 




Earth 
Asteroid 
Trajectory 
SOI Earth 


(b) 2 nd through 4 th Stages, Sun-Centered 


Figure 3. Trajectory Schematics 


We designed the algorithm to use an embedded trajectory grid methodology to assess NEA ac- 
cessibility by computing all the possible round-trip trajectories. A schematic of the parameter space 
structure of the algorithm is depicted in Figure 4. The algorithm varies six independent parameters 
within set bounds in order to identify solutions that satisfy problem constraints. For this paper, we 
chose to use the same constraints as in Reference 4 so that results could be compared. Whether an 
NEA is classified as accessible (from a dynamics point of view) depends on maximum allowable 
round-trip flight time, lower and upper bounds for the trajectory segment flight times, step sizes at 
which trajectory segment flight times are sampled, lower and upper bounds for the Earth departure 
date, the step size at which Earth departure date is sampled, and constraints on total spacecraft mass 
as a function of the assumed launch and propulsion system capabilities. The independent parameters 
are: 

• tdepe : D ate of Earth Departure 

• #o: The initial angle at LEO departure 

• a to : The location of the Moon at departure time, in accordance with Eq. (6) 


7 


• t°farr a : Time of flight from SOI e to Asteroid 

• to f stCL y : Time to stay at NEA 

• t°f arVe : Time of flight from Asteroid to Earth 



Figure 4. Schematic of Embedded Grid Algorithm 

To evaluate whether the total mission duration constraint is met, the total mission duration is 
computed as the sum of the flight times associated with each of the four mission stages, according 
to 

to f tot — to f soi + tof arra + to f s tay H - tof arTe < 360 days. (13) 

Note that tof so i, the time of flight to get from LEO to SOI e depends on the specific values chosen 
for #o and a to . 

Evaluation of mission mass constraints requires the computation of the total mission AV, which 
is simply the sum of all the maneuver magnitudes, 

AVtot — AV S oi + AV arra + A Vdep a T AV arTe . (14) 

In this study, we utilized the launch mass versus C 3 values of the notional Ares V heavy-lift launch 
vehicle. We assumed that the capabilities of Ares V is representative of a future heavy-lift human- 
rated launch vehicle that will be available by the mid-2020s . 7 For a C 3 = —2.0390 km 2 /s 2 , this 
corresponds to an available launch mass mc 3 = 5.5774 • 10 4 kg. The mass constraint is satisfied if 

^ < 1 (15) 

mc 3 
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where 


( AV tot \ 

Wlreq — J ^ 

and AVtot is in (14). The parameters nidry = 17,078 kg and Isp = 314 s are obtained from 
the notional parameters of the Orion crew vehicle, which was chosen because it is the most recent 
design capable of carrying crew on a mission beyond LEO . 4 

Preliminary Results: NEA 2000 SG 344 

NEA 2000 SG 344 has an absolute magnitude, H , of 24.8, which yields an estimated physical 
size of 29 to 66 m for an albedo range of 0.25 to 0.05, respectively. While those albedos are 
typical for NEAs, it is quite possible for an NEA to have a significantly lower or higher albedo 
and thus the estimated size based on absolute magnitude and assumed albedo should be treated as 
very approximate. 2000 SG 344 has semi-major axis a = 0.977 AU, eccentricity e = 0.066, and 
inclination i = 0.11°. With such an Earth-like orbit, one would expect this NEA to offer many 
efficient round-trip trajectory opportunities when in proper phase with respect to Earth. 

Indeed, the direct approach in Reference 4 identified over 700, 000 trajectory solutions that met 
the given constraints on Earth departure date, total round-trip flight time, and mission mass budget. 
In Reference 4 the non-dimensional parameter called the mass ratio, denoted by a m , was developed 
as a criterion for accessibility. a m , defined in Equation (15), is equal to the required initial total 
spacecraft mass (dry mass plus propellant) divided by the payload mass capability of a notional Ares 
V heavy-lift vehicle. The required propellant load for the spacecraft is a function of the required total 
mission AV for a particular trajectory solution while the payload mass capability of the notional 
Ares Visa function of the Earth departure C 3 for the particular trajectory solution. In Reference 4, 
an NEA had to offer at least one trajectory solution (within given constraints) for which a m < 1; 
2000 SG 344 offered over 700, 000 such solutions, with Earth departure years between 2027 and 
2030. 

The new algorithm described herein was exercised on 2000 SG 344 and produced nearly 4, 400 
trajectory solutions that include a lunar flyby in the outbound trajectory segment departing Earth. 
However, these solutions occur on only four unique departure dates that fall between the years 2028 
and 2030. It is important to note that only solutions which required AV SO i < 0.01 km/s were 
considered viable. This was to ensure that the outbound trajectory can be essentially smooth, with 
no major maneuver requirements. Table 1 compares the results of the new algorithm to the original 
results in Reference 4. Note that the restriction of including a lunar flyby dramatically reduced 
the number of viable trajectory solutions, reduced the window of viable of Earth departure dates 
from three years to two years, and increased the minimum available total round-trip flight time from 
three months to four months. However, also note that the minimum available mass ratio for the lunar 
flyby cases is about 2% lower than for the direct method. This provides our first demonstration that 
the inclusion of a lunar flyby in the outbound trajectory segment can indeed reduce the required 
propellant mass for the mission. 

An example trajectory for 2000 SG 344 with initial condition given in Table 2 is shown in Fig- 
ures 5(a) and 5(b) in a Sun-centered frame and in an Earth centered frame in Figures 5(c) and 5(d). 
Table 2 compares this specific trajectory solution generated by the new lunar flyby method to one 
generated by the direct method of Reference 4. Note that these trajectory solutions have nearly iden- 
tical Earth departure dates (one day apart) and total round-trip flight times (151 versus 154 days). 
However, for the lunar flyby solution, the Earth departure C 3 is considerably lower (and is in fact 
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Table 1. Comparison of Lunar Flyby and Direct Trajectory Scan Results for NEA 2000 SG 344 



Lunar Flyby 

Direct 

# of Viable Trajectory Solutions 

4,398 

708, 703 

Departure date 

2028 - 2030 

2027 - 2030 

f min (days) 

128 

96 

f max (days) 

360 

360 

Min mass ratio a m 

0.3429 

0.3630 

Max mass ratio a m 

0.9994 

1.0000 


commensurate with typical C 3 requirements for lunar missions), the required AV at Earth’s Sphere 
of Influence is extremely small, and the total mission AV requirement is significantly reduced. Thus 
the mass ratio for the lunar flyby case is lower since the available launch vehicle payload capacity 
is higher at the lower C 3 , and the total required spacecraft propellant mass is less due to the reduced 
total mission AV (all else being equal). 

Table 2. Comparison of Lunar Flyby and Direct Trajectory Solutions for NEA 2000 SG 344 



flyby 

direct 

8o n 

0.00 

- 

Q to (°) 

1 15.00 

- 

Launch date 

01/09/2029 

01/10/2029 

tof (days) 

151 

154 

tofsoi (days) 

11 

- 

tofarr a (days) 

68 

76 

tof stay (days) 

20 

20 

tofarre (days) 

52 

58 

C 3 (km 2 /s 2 ) 

-2.0390 

1.6579 

A V soi (km/s) 

0.0068 

- 

AVtot (km/s) 

3.4312 

4.5432 

mass ratio a m 

0.9331 

0.9358 


While these results are certainly encouraging, recall that thus far two major assumptions have 
been made. First, the trajectory of the Moon is modeled as being circular about the Earth, and in 
the ecliptic plane. Second, the lunar flyby is constructed with the freedom to place the Moon at 
precisely the correct location at the necessary time. However, subsequent analysis of the true lunar 
ephemeris revealed that Moon is never at the proper location on any of the four available departure 
dates. An example of this situation is shown in Figure 6 and clearly demonstrates that the Moon is 
rather far from the required location. Therefore, a round-trip trajectory utilizing a lunar flyby for 
Asteroid 2000SG344 is not possible in actuality. However, the results obtained with the fictitious 
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Trajectory in Sun Centered Frame 


x iq 8 Trajectory in Sun Centered Frame 

Reach SOI - 


O Sun 

Earth 

- Asteroid 
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(a) Sun-Centered, Ecliptic Plane (b) Sun-Centered, Inclined. Axis are Not at Same Scale 



(c) Earth-Centered, Lunar Flyby Close-Up 


x iq 6 Trajectory in Earth Centered Frame 



(d) Earth-Centered, Round-Trip Trajectory 


Figure 5. Example Theoretical Trajectory for NEA 2000 SG 344 . 


lunar ephemeris served as a strong proof of concept and motivated the continuation of this research 
using the true lunar ephemeris. 

MOON’S EPHEMERIS IMPLEMENTATION 

The true lunar ephemeris was incorporated into our lunar flyby algorithm following the successful 
proof of concept for NEA 2000 SG 344 using the fictitious lunar ephemeris. The previous methodol- 
ogy utilizing the fictitious model for the Moon’s orbit still serves to provide good initial guesses for 
outbound trajectories that be may be able to utilize a lunar flyby en-route to a NEA. Thus the first 
step in the upgraded algorithm still involves constructing the trajectory sequence using the fictitious 
Moon and then checking the Moon’s actual ephemeris to see whether its true position is within some 
tolerance of where the theoretical model indicates it needs to be on the particular flyby date. If so, 
then the initial guess can be used, in combination with the true lunar ephemeris, to compute the true 
trajectory solutions. 

It is worth noting that the Moon’s true state, x, is in R 3 , which means that a lunar flyby can 
provide some inclination change for the spacecraft; modeling this was not possible with the planar 
model for the Moon’s orbit. Also, the actual orientation of the Moon’s orbit plane with respect to 
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Figure 6. Moon’s location is not coincident with trajectory. 


the ecliptic (a difference of about 5°) is now accounted for when transforming coordinates from the 
Earth-centered to the Sun-centered frame. 

The Case of NEA 2003 LN 6 

Since NEA 2000 SG 344 was previously determined to have no lunar flyby opportunities when the 
true lunar ephemeris is accounted for (and since it was already accessible in the previous study with- 
out a lunar flyby), we selected a different candidate target from the list of marginally inaccessible 
NEAs presented in Reference 4. The goal was to show that a lunar flyby alone could make the NEA 
accessible according to the constraints and criteria of the previous study. We selected NEA 2003 
LN 6 in the hopes that the true lunar position would be favorable for lunar flyby trajectory solutions. 
This NEA has an estimated size of 35 to 77 m, its orbit has a semi-major axis of 0.857 AU, its 
orbital eccentricity is 0.211, and its orbital inclination is 0.633°. The best-case (minimum) mass 
ratio trajectory solution found in the previous study for this NEA had an associated mass ratio value 
of 1.00325. Thus the required spacecraft propellant mass was just beyond the payload capacity of 
the assumed launch vehicle. 

The first step was to employ the algorithm as described previously, using the fictitious model for 
the Moon’s orbit. This produced 36 possible round-trip trajectory solutions which included a lunar 
flyby. These possible trajectory solutions were spread over four departure dates within the year 
2025. Again, note that only solutions with AV^ < 0.1 km/s were considered. Table 3 compares 
the results from the initial lunar flyby trajectory scan to the results of the direct method trajectory 
scan given for the previous study in Reference 4. Note that in the previous study only one minimum 
mass ratio trajectory solution was tabulated for each NEA found to be inaccessible. Note also that 
even using the fictitious model for the Moon’s orbit, we already reduce the spacecraft propellant 
mass sufficiently to obtain an a m < 1 solution for this NEA. 

Out of the 36 possible trajectory solutions, 8 were found to have the Moon actually near the 
required location. One of these trajectory solutions was selected and the true lunar ephemeris was 
then used to compute a final lunar flyby trajectory solution for NEA 2003 LN 6 . The details of 
this trajectory solution are presented in Table 4. Note that the restriction of C 3 in Eq. ( 8 ) having 
to be that of a Hohmann transfer is absent in the model with the Moon’s true ephemerides, since 
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Table 3. Comparison of Lunar Flyby and Direct Trajectory Scan Results for NEA 2003 LN6 



Lunar Flyby 

Direct 

# of Viable Trajectory Solutions 

36 

0 

Departure date 

2025 

2025 

TOFmin (days) 

188 

204 

TOFjnax (days) 

196 

- 

Min mass ratio 

0.9929 

1.00325 


C 3 = -2.0252 km 2 /s 2 as shown in Table 4. The associated trajectory plots are presented in 
Figure 7. The comparison between the theoretical and true trajectories is shown in Figure 7(e). 

Table 4. Lunar Flyby Solution for NEA 2003 LN6 Utilizing True Ephemerides Throughout 


Quantity 

Value 

Est. Size (m) 

35-77 

Launch date 

12/09/2025 

#0 

331.61° 

C z (km 2 /s 2 ) 

-2.0252 

TOF tot (days) 

198.4 

tofsoi (days) 

11.91 

tofdep (days) 

134.49 

to f stay (days) 

4 

to fret (days) 

48 

AV soi (km/s) 

0.1122 

A V t ot (km/s) 

3.6428 

Q m 

0.9997 
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(a) Sun Centered, Ecliptic Plane 



x io 6 Trajectory in Sun Centered Frame 



(b) Sun Centered, Inclined. Axis are Not at Same Scale 


x io 6 Trajectory in Earth Centered Frame 



(c) Earth Centered, Flyby 


(d) Earth Centered 


Trajectory in Earth Centered Frame 



(e) Comparison of Trajectory Solutions Using Theoretical and True Lunar 
Ephemerides 

Figure 7. Plots of Round-trip Trajectory to 2003 LN6 Using True Lunar Ephemeris 
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CONCLUSION 


In this paper we present a new general methodology for designing round-trip trajectories to NEAs 
in which lunar gravity assist is utilized to aid the outbound trajectory segment. Our results were 
compared to those from a previous study 4 in which round-trip trajectories were designed without 
the lunar gravity assist. The objective was to gain sufficient velocity from Moon’s gravity to reduce 
the total spacecraft propellant mass required for a round-trip mission to an NEA. Asteroid 2003 LN 6 , 
an NEA categorized as marginally inaccessible in Reference 4, was made accessible (according to 
the previous study criteria) via a lunar flyby using the new algorithm presented herein. The lunar 
flyby method accomplished this by reducing the post-Earth departure AV for the mission by only « 
2%, while also reducing the Earth departure C 3 to be commensurate with what is typically required 
for a lunar mission. 

The new algorithm minimizes the magnitude of the post-lunar flyby maneuver required to reach 
the NEA, AVgoi- In future work we will seek to eliminate the need for that maneuver, i.e., achieve 
AV so i = 0. We also plan to optimize the trajectories such that that both AV tot and iof tot are mini- 
mized, while also converting all the impulsively modeled maneuvers to finite bums. The possibility 
of a lunar flyby in the return segment from an NEA of the trajectory will be studied as well to 
determine if other NEAs can be made accessible that are now classified as inaccessible. 
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